Anomalous spin Hall effects in Dresselhaus (110) quantum wells 
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Anomalous spin Hall effects that belong to the intrinsic type in Dresselhaus (110) quantum wells are dis- 
cussed. For the out-of-plane spin component, antisymmetric current-induced spin polarization induces opposite 
spin Hall accumulation, even though there is no spin-orbit force due to Dresselhaus (110) coupling. A surpris- 
ing feature of this spin Hall induction is that the spin accumulation sign does not change upon bias reversal. 
Contribution to the spin Hall accumulation from the spin Hall induction and the spin deviation due to intrinsic 
spin-orbit force as well as extrinsic spin scattering, can be straightforwardly distinguished simply by reversing 
the bias. For the inplane component, inclusion of a weak Rashba coupling leads to a new type of S y intrinsic 
spin Hall effect solely due to spin-orbit-force-driven spin separation. 

PACS numbers: 72.25.Dc, 71.70.Ej, 73.23.Ad 



I. INTRODUCTION 

The intensive efforts on spin Hall effect (SHE) both exper- 
imentally and theoretically during the past decade have suc- 
cessfully built another milestone in condensed matter physics. 
Spin separation in semiconductors is not only possible but 
natural, so that manipulating spin properties of charge car- 
riers in electronics is promising. The earliest theoretical idea 
that up and down spins may laterally separate upon transport 
due to asymmetric scattering was proposed in 1971. 1 - 2 More 
than three decades later, the power of optical measurements 
on high quality mesoscopic samples made SHE in semicon- 
ductors no longer an idea but an experimental fact. 3 Right 
before the first observation of Ref. 3 in 2004, mechanisms 
of SHE was further extended from spin-dependent scattering 
that was later categorized as extrinsic, to spin-orbit-coupled 
band structure that was later categorized as intrinsic. 4 ' 5 Ex- 
perimentally, most observations so far have been attributed to 
the extrinsic SHE, 3 ' 68 while evidence of the intrinsic SHE 9 is 
relatively few. Nonetheless, intrinsic SHE remains an impor- 
tant issue that until now still receives enduring efforts. 10 

In the intrinsic SHE, spin separation is solely due to the 
underlying spin-orbit coupling in the band structure, so that 
SHE can exist even in systems free of scattering (but of finite 
sizes 11,12 ). In the ballistic limit, the spin separation can be 
vividly visualized by the transverse spin-orbit force 13 ' 14 de- 
rived by using the Heisenberg equation of motion, 
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where r is the position operator and T-L is the single- 
particle Hamiltonian. For well discussed two-dimensional 
systems with Rashba coupling 15 described by Hr = 
(a/h) (pya x — p x o~ v ), as well as linear Dresselhaus (001) 



coupling 16 ' 17 described by Tip 1 
the spin-orbit force is given by 13 
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pling constants, respectively, and <j x 1 a v 1 a z are Pauli matri- 
ces. Equation (2) clearly depicts a lateral spin deviation of the 
S z = (h/2)a z spin component with opposite contributions 
from Rashba and Dresselhaus (001) couplings, as sketched in 
Fig. 1(a) and (b), respectively. 

SHE in Dresselhaus (110) quantum wells (QWs), on the 
other hand, is relatively less discussed theoretically, 18 al- 
though a previous experimental effort 6 had revealed in GaAs 
(1 10) QWs the existence of SHE that was attributed to the ex- 
trinsic type. In this paper, anomalous SHEs that belong to the 
intrinsic type in Dresselhaus (110) QWs are discussed. We 
show that the spin Hall pattern of S z can be induced when 
the transport direction is properly oriented, even though the 
Dresselhaus (110) coupling does not result in spin-orbit force 
to separate opposite S z spins upon transport [see Fig. 1(c)]. 
Moreover, we propose a Rashba-coupling-assisted intrinsic 
SHE in S y that is truly due to spin-orbit force under the in- 
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where p is the momentum and e z is the unit vector of the 
plane normal. Here a and /3 are Rashba and Dresselhaus cou- 



FIG. 1. (Color online) Spin deviation in S z due to spin-orbit force 
in (a) Rashba and (b) linear Dresselhaus [001] systems, (c) Spin 
Hall induction due to antisymmetric CISP along [001] in Dressel- 
haus [110] systems may cause opposite S z accumulations as well, 
even though there is no spin-orbit force, (d) In the case of linear 
Dresselhaus [110] plus weak Rashba couplings, a coupled spin-orbit 
force given by Eq. (11) may lead to spin deviation in S y . 



teraction of Dresselhaus (110) plus a weak Rashba couplings 
[Fig. 1(d)]. 

This paper is organized as follows. In Sec. II we briefly 
introduce the formulas required in the Landauer-Keldysh for- 
malism employed in the numerical analysis of Sec. Ill, where 
we visualize the proposed spin Hall induction and Rashba- 
coupling-assisted SHE in S y . Comparison of the present bal- 
listic calculation with the diffusive experiment of Ref. 6 will 
be discussed and the transport parameters used in our numer- 
ical data will be remarked. We conclude in Sec. IV. 



II. FORMULAS 
A. Linear Dresselhaus (110) coupling 

The Dresselhaus (110) coupling up to the term linear in mo- 
mentum can be written as 
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where x, y, and z axes are chosen along [110], [001], and 
[110], respectively. Throughout the present discussion, we 
will focus on this Dresselhaus (110) linear term, so that 
without ambiguity we use the same notation f3 to denote its 
coupling strength. The spin-orbit field subject to Eq. (3) 
is depicted in Fig. 2(a). Clearly, when propagating with 
k = (±\k x \ ,k y ) electrons encounter antisymmetric spin- 
orbit fields on the left and right sides of the [001] axis (or the 
y axis). Hence the current-induced spin polarization (CISP) 
effect 19-21 is expected to build opposite S z spin densities at 
the two sides, as conceptually depicted in Fig. 1(c). 

To better illustrate this spin Hall induction, we will in 
Sec. Ill first consider a T-bar ballistic nanostructure, attached 
to left, right, and bottom leads (from the top view) that are 
made of normal metals. The central region is described by the 
square-lattice tight-binding Hamiltonian, 
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where the sum over (nm) of the second term is run for the 
sites nearest to each other, satisfying |r m — r„| = a, a being 
the lattice grid spacing and r n the position vector of site n, 
and the hopping matrix is given by 



= -£ 1 - it D d x a z 
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Here U is the on-site energy set to be constant over the whole 
sample, £ = h 2 /2ma 2 is the kinetic hopping energy, 1 is 
the identity, c m (c^) annihilates (creates) an electron at site 
m, tf) = f3/2a is the Dresselhaus hopping parameter, and 
d x = (r, n — r n ) • e x is the hopping displacement along x 
from site m to site n. 



B. Landauer-Keldysh formalism 

To image the nonequilibrium charge, charge current, spin, 
and spin current densities under the influence of the biased 



FIG. 2. (Color online) (a) Linear Dresselhaus [110] spin-orbit field 
in k space, (b) Local charge density (eJV) (color shading) and local 
charge current density (arrows) in a 16a x 8a Dresselhaus [1 10] sam- 
ple subject to three terminals with a weak bias eVo = 10 _3 to- (c) 
Local spin density (<j z ) (color shading) and its corresponding local 
spin current density (J ) (arrows) in the same device as (b) with 
same conditions. Spin Hall induction in a 40a x 40a Dresselhaus 
[110] sample with (d) upward bias and (e) downward bias; strong 
bias eVo — 0.4io is applied. Note that in (b)-(e), the regions outside 
the dashed lines are simulating the leads (zero spin-orbit coupling 
and constant on-site energy set equal to the applied bias). 



leads, a powerful and convenient approach is the Landauer- 
Keldysh formalism, 22,23 especially for the present ballistic 
case free of particle-particle interaction. In this formalism, 
physical quantities in a nonequilibrium but steady state are 
expressed in terms of the lesser Green function matrix G < , 
provided that those physical observables of interest are well 
defined. 24 Each matrix element G,f m in our spin- 1/2 electron 
system is a 2 x 2 submatrix, so that the size of full G < amounts 
to 2N x 2N, N being the total number of lattice grid points. 
Following Ref. 24 with moderate extension, we have and will 
use the local charge and spin densities 
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for the flow from site n to site m. Here Tr s is the trace done in 
the spin space, the explicit energy E dependence of G< n (E) 



is suppressed, and {A, B} = AB + BA is the anticommuta- (a) 

tor. For the present illustration of the spin Hall induction, we 
will first consider a pure linear Dresselhaus (110) system and 
use the hopping matrix (5) with to = O.Uq, which is within a 
reasonable range. Parameters extracted from experiments will 
be discussed later in Sec. HID. Other transport parameters are 
as follows: hopping parameter to = 1, on-site energy U = 
(so that band bottom Ef, = 0), Fermi energy is 0.2£o above 
Et (so that the square lattice simply serves as the grid of a 
free electron gas). We will always label + and — to indicate im 

an applied bias voltage of +eVb/2 and — eVo/2 on each lead, 
respectively, with eVb > 0. (Note that e = — |e| is the neg- 
ative electron charge, and hence electrons always flow from 
+ to — signs). In the rest of our analysis, we will focus on 
the nonequilibrium contribution 24 of those quantities listed in 
Eqs. (6)-(9), and hence the integration range will be taken as 
E F - eV a /2 ->• E F + eV /2. 



FIG. 4(a) 
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III. NUMERICAL ANALYSIS 
A. Spin Hall induction in S z : Pure Dresselhaus (110) coupling 

Employing the Landauer-Keldysh formalism briefly intro- 
duced above, we now drive electrons in the T-bar nanostruc- 
ture from bottom to left and right leads with eVo = 10 _3 io, 
as shown in Fig. 2(b), where the background color shading is 
determined by the local charge density (eN) given in Eq. (6), 
while each arrow indicates the local charge current density 
given by Eq. (8). In Fig. 2(c), the color shading is determined 
by the local spin (S z ) density [Eq. (7)] and clearly shows an 
antisymmetric S z polarization for electrons moving into the 
left and right leads because of the opposite Dresselhaus (110) 
fields they feel. The local spin current density indicated by the 
arrows therein is given by Eq. (9), which is derived from the 
symmetrized spin current operator J 5, = {J m ^ n /e,Si} /2 
in a way similar to Ref. 24. A pure spin current from right 
to left is observed; at right side the spin current is flowing to- 
ward left because of the negative S z times the right moving 
particle current while at left side the left flowing spin current 
stems from the product of the positive S z and the left moving 
particle current. 

For a two-terminal device made of Dresselhaus (110) QW 
oriented along [001] (y axis), spin Hall accumulation of op- 
posite 5*2 is therefore expected, as shown in Fig. 2(d)-(e), 
where we consider a strong bias voltage of eVo = 0.4io for 
a 40a x 40a sample. A striking difference between the spin 
Hall induction introduced here and the spin Hall deviation due 
to spin-orbit force is the independence of the accumulation 
sign on the bias direction. Whether driving the electrons from 
bottom to top [Fig. 2(d)] or from top to bottom [Fig. 2(e)], 
one always observe a negative S z accumulation at right while 
positive at left. 

Contrary to the present ballistic nanostructure here, the 
SHE previously observed in GaAs (110) QWs 6 was in dif- 
fusive regime and attributed to the extrinsic type. The experi- 
ment used ac lock-in detection referenced to the frequency of 
a square wave alternating voltage with zero dc bias offset, and 



FIG. 3. (Color online) Imaging of local spin densities (a z ) and (o v ) 
in a 60a x 60a [110] sample with (a)-(b) top- to-bottom and (c)- 
(d) right- to-left orientations. A strong Dresselhaus [110] coupling 
to — 0.08£o and a weak Rashba coupling tn — 0.02£o are used. 
Bias is set eVb = 0.4io- 



the resulting signal is sensitive only to the difference in Kerr 
rotation between positive and negative bias. Our spin Hall in- 
duction that does not depend on bias direction may therefore 
either hardly contribute or be subtracted in the result of Ref. 
6. 



B. Spin Hall deviation in S y : Strong Dresselhaus (110) with 
weak Rashba couplings 

Next we recall the spin-orbit force Eq. (1). For pure Dres- 
selhaus (110) systems given by Eq. (3), there is no way to 
obtain a nonvanishing F so since eventually the Pauli matrix 
a z will commute with itself, even if the cubic term that is still 
in terms of a z is involved. The only possibility in this case for 
a nonvanishing F so to survive is to introduce spin-orbit terms 
involving a x or a v . Combination of Rashba coupling with 
the present linear Dresselhaus (110) term is therefore a natu- 
ral candidate, which is possible for, for example, asymmetric 
GaAs (110) QWs, as are the cases of Ref. 6. The spin-orbit 
force for this Rashba-Dresselhaus (110) QW is 
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Without Rashba term a, the spin-orbit force vanishes and zero 
spin current is hence expected. From a gauge viewpoint, the 
existence of equilibrium spin current in (110) QWs will re- 
quire Rashba term to break the pure gauge. 25 Note also that 
the a squared dependence for the a z component in Eq. (10) is 
similar to the result in Ref. 26. 

Here of particular interest is the case of weak Rashba cou- 
pling, such that Eq. (10) becomes 
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which predicts a lateral spin Hall deviation in S y that re- 
quires a weak but nonzero Rashba coupling a. To further 
visualize the predicted Rashba-assisted S y SHE, we con- 
sider a 60a x 60a sample with Dresselhaus (110) hopping 
to = 0.08io and Rashba hopping tn = a/2a = 0.02io, 
attached to two leads under a bias voltage eVo — 0.4£rj. For 
the [001] -oriented (electron flow along —y) sample, the S z 
spin Hall pattern due to spin Hall induction is observed in Fig. 
3(a). Meanwhile, an S y spin Hall pattern is also shown in 
Fig. 3(b), which is a combined consequence of not only the 
spin deviation Eq. (11) but also an antisymmetric CISP by the 
Rashba coupling. Along the — y axis, electrons with wave vec- 
tor k = (± \k x \ , — \k y \) encounter opposite y component of 
the clockwise Rashba spin-orbit field: negative for + \k x \ and 
positive for — \k x \. Hence a spin Hall induction in S y due to 
Rashba coupling contributes to Fig. 3(b) as well. In addition, 
the contribution of the spin-orbit force Eq. (1 1) predicts a +S y 
(—S y ) accumulation at left (right) side of the electron flow, 
for lateral distance shorter than the spin precession length L so 
(around 15a here); the accumulation sign reverses when the 
lateral distance exceeds L so , as is the case in our 60a x 60a 
here. Therefore the two contributions, spin Hall induction and 
spin-orbit force Eq. (11), are additive in Fig. 3(b). 

For the [110]-oriented (electron flow along — x) sample, 
there is a vague spin Hall pattern in S z because of the ab- 
sence of the antisymmetric CISP and weak spin-orbit force 
[Fig. 3(c)]. The average of (S z ) over the whole sample ba- 
sically reveals the usual CISP effect as observed in Ref. 6. 
The S y pattern, on the other hand, exhibits a clear spin Hall 
accumulation pattern which is solely attributed to the spin- 
orbit force Eq. (11), as shown in Fig. 3(d). As explained, — S y 
(+S y ) accumulates at left (right) side of the electron flow be- 
cause the lateral distance has exceeded L so . Upon the bias 
reversal, the ±S y edge accumulations swap (not shown here), 
which is the general feature of the spin Hall pattern due to spin 
deviation by intrinsic spin-orbit force, as well as by extrinsic 
spin scattering. The spin Hall induction such as that of S z in 
the Dresselhaus (110) case along ±y, however, does not have 
this feature. Another difference between the spin Hall patterns 
induced by antisymmetric CISP and spin-orbit-driven spin de- 
viation is that in the former the signs of the spin accumulation 
do not change with the increasing sample width, while in the 
latter they do. This difference can also be told in Fig. 3: con- 
stant sign in each lateral side in panel (a) but varying sign in 
panels (b) and (d). 



C. From pure Dresselhaus (110) to pure Rashba cases 

Finally, we laterally scan the local spin densities (S z ) and 
(S y ) in Fig. 4 at the positions marked by the dashed lines in 
Fig. 3, for a set of various spin-orbit coupling parameters from 
pure Dresselhaus (110) (black curves) to pure Rashba (light- 
est gray curves). For the [001] -oriented sample, S z spin Hall 
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FIG. 4. (Color online) Local spin densities (<r z ) and {a v ) in a 
60a x 60a sample as a function of the transverse position with var- 
ious coupling parameters to and tR. The origin is set at the cen- 
ter of the sample. As indicated in panel (a), the coupling param- 
eters from black to the lightest curves are (tR,to) = (0,0.1) to, 
(0.02, 0.08) to, (0.04, 0.06) to, (0.06, 0.04) to, (0.08, 0.02) to, and 
(0.1, 0) to- In each panel, the red (dark gray) thick curve correspond 
to Fig. 3. 



pattern shown in Fig. 4(a) gradually evolves from spin Hall in- 
duction due to Dresselhaus (110) coupling to spin Hall devia- 
tion driven by spin-orbit force due to Rashba coupling. In Fig. 
4(b), the black curve for the pure Dresselhaus (110) shows 
zero everywhere [and so are those for Fig. 4(c)-(d)], while 
a weak Rashba coupling assists the formation of the S y spin 
Hall pattern; the antisymmetric pattern holds all the way to 
pure Rashba because the Rashba coupling contributes through 
the antisymmetric CISP in this orientation as explained previ- 
ously. For the [110]-oriented sample, turning on of the weak 
Rashba coupling builds S y spin Hall pattern [Fig. 4(d)] but 
not too much for S z [Fig. 4(c)]. Down to pure Rashba, the S z 
pattern recovers the spin Hall accumulation due to spin-orbit 
force [Fig. 4(c)], while that for S y shows symmetric CISP 
[Fig. 4(d)]. 



D. Remark on transport parameters 

In our numerical analysis for the pure Dresselhaus case, we 
have set to /to = 0.1 mostly based on an illustrative reason. 
This coupling ratio allows a direction comparison with Ref. 
22, where tn/to — 0.1 is chosen, in the later part of coexisting 
Rashba and linear Dresselhaus (110) couplings (such as Fig. 
4). 

Comparing with the GaAs (110) QWs used in the exper- 
iment of Ref. 6, the coupling ratio to /to mav t> e one order 
weaker than ours. The sample they used behaves like a single 
75 A Al .iGa .gAs QW. Using the relation /3 = l{k 2 z ) with 
hard wall approximation (k 2 ) w (ir/w) and 7 « 27eVA 
for both GaAs and InAs QWs, 27 this well width of w = 75 A 



TABLE I. 

Ref. 27. 
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mass and Dresselhaus coefficients taken from 


QW type 


GaAs 


AlAs InAs InSb CdTe ZnSe 


m/mn 
7 (eVA 3 ) 


0.067 
27.58 


0.15 0.023 0.014 0.09 0.16 
18.53 27.18 760.1 43.88 14.29 



In general for a stronger to /to, which can be rewritten as 
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FIG. 5. (Online color) Coupling ratios to /to estimated by Eq. (12) 
for various QWs with widths w — 5 nm, 7.5 nm, 10 nm. 



leads to/3w 4.74 x 10 -2 cV A. Effective mass was reported 
to be m = 0.074too, ttiq the electron rest mass. The sheet 
density is n s = 1.9 x 10 12 cm -2 , which allows us to estimate 
the location of the Fermi energy 23 Ep — Ef, = Trh 2 n s /iri « 6. 
15 x 10~ 2 eV. In order for the long wavelength limit to 
be valid, the chosen lattice constant a has to yield a kinetic 
hopping constant to that keeps Ep close to E^. Choosing 
a = 2 nm leads to t » 0.13 cV so that Ep - E b « 0.l2t is 
satisfying (recall Ep —Ef, = 0.2to in our numerical results as 
well as in Refs. 22 and 24). The coupling ratio with this a is 
to/to ~ 0.01. The Rashba strength in Ref. 6 was reported to 
bea = 0.018 cV A, leading to a//3 « 0.38, which is not too 
far from our t^/tp, = 0.25 in Sec. IIIB. Replacing these pa- 
rameters in our results does not change significantly the main 
features we have shown. 

Coupling ratio of tp>/to — 0.1 is actually possible for 
QWs with stronger Dresselhaus bulk coefficient 7. For InSb 

QWs, 27 we have 7 = 760 eV A . Consider the Inn.sgGao.nSb 
QWs with effective mass m = 0.018r7in and sheet electron 
concentration n s = 2.9 x 10 11 cm" 2 reported in Ref. 28, 
where the QW width is relatively thick: 30 nm. If reducing 
the QW width to 7.5 nm, which is common in GaAs QWs, 
and assuming a = 3 nm, the coupling ratio is estimated as 
to/to ~ 9. 45 x 10" 2 , close to our to/to = 0.1. The Fermi 
energy in units of to is (Ep — Eb)/to = 2na 2 n s w 0.16, 
which is also close to our (Ep — Eb)/to = 0.2. Hence the 
transport parameters used in our calculation are within a rea- 
sonable range. 
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a larger product 7717, and a thinner QW width w will be re- 
quired. The effective mass m and Dresselhaus coefficient 
7 for various QWs taken from Ref. 27 are collected in Ta- 
ble I, For these QWs we use Eq. (12) with a = 3nm to 
summarize the coupling ratio tp>/to in Fig. 5 for QW widths 
w = 5 nm, 7.5 nm, 10 nm. 



IV. CONCLUSION 

In conclusion, we have shown that spin Hall induction for 
±[001] transport in (110) QWs due to antisymmetric CISP of 
linear Dresselhaus coupling that yields zero spin-orbit force is 
possible to generate a spin Hall accumulation pattern in S z , 
whose signs do not depend on the bias direction. Experimen- 
tal investigations with a dc bias offset in ballistic (or at least 
quasi-ballistic) III-V (110) symmetric QWs may potentially 
identify our proposed effect. From the coupling ratios summa- 
rized in Fig. 5, InSb (110) QW is promising for the presently 
proposed spin Hall induction in S z , while InAs is less sug- 
gested. A new type of spin Hall deviation in S y is also pre- 
dicted in the Dresselhaus (110) QWs in the presence of a weak 
Rashba coupling. Experimental observation for this S y spin 
Hall effect may require a good control over the Rashba and 
Dresselhaus couplings, which has been proved possible for 
(001) QWs, 2931 and should be achievable also for (110) QWs. 
We categorize these two intrinsic spin Hall mechanisms — spin 
Hall induction in S z and spin Hall deviation in S v , as anoma- 
lous SHEs. 

We note that the Dresselhaus cubic term, neglected in the 
present study, will become important when Fermi wave vec- 
tor kp is long or QW width w is thick. In this case the spin 
Hall induction in S z along ±[001] axis as discussed above 
should remain, while additional spin Hall induction axes close 
to [111] and [112] will further emerge; see Fig. 6.20 in Ref. 
27 and Fig. 2(a) in Ref. 6. Inclusion of the Dresselhaus cubic 
term is left as a future extending work. 
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